function [ b ] = assemblyBoundaryVector( domain2D )
%ASSEMBLYBOUNDARYVECTOR Summary of this function goes here
%   Detailed explanation goes here

n = size(domain2D.e1, 1);

b = zeros(2*n, 1);

for i = 1 : n
   i1 = domain2D.e1(i);
   i2 = domain2D.e2(i);
   
%    if((domain2D.vx(i1)>-0.54)&(domain2D.vx(i2)>-0.54)&(domain2D.vx(i1)<-0.5)&(domain2D.vx(i2)<-0.5)&(domain2D.vy(i1)>0.9)&(domain2D.vy(i2)>0.9))
%       b(2*i - 1) = -10;
%    end
%    
%    if((domain2D.vx(i1)>-0.3)&(domain2D.vx(i2)>-0.3)&(domain2D.vx(i1)<-0.25)&(domain2D.vx(i2)<-0.25)&(domain2D.vy(i1)>0.9)&(domain2D.vy(i2)>0.9))
%       b(2*i - 1) = 10;
%    end

   if((domain2D.vx(i1)<-0.999)&(domain2D.vx(i2)<-0.999))
      b(2*i - 1) = -5;
   end
   
   if((domain2D.vx(i1)>0.999)&(domain2D.vx(i2)>0.999))
      b(2*i - 1) = 5;
   end

end




% br = reshape(b', n, 2);

% br((domain2D.vx>-0.01)&(domain2D.vx<0.01)&(domain2D.vy>0.9),2) = 1;

% b = reshape(br, 2*n, 1);

% b(6) = 1;

% b(237) = 5;

end

